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1. Abstract 

We study the stochastic quantization of two-dimensional nonlinear sigma model in 
the large N limit. Our main tool is the effective Langevin equation with which we 
investigate nonperturbative phenomena and derive the results which are same as the 
path integral approach gives. 



Recently, the stochastic quantization method' 11 [21 has been used to study sponta- 
neous (dynamical) breakdown of symmetry,' 31 for example, Nambu-Jona-Lasinio model.' 41 Their| 
main tool for investigation is the stochastic generating functional from which the effec- 
tive potential is derived. On the contrary, the Langevin equation, which is very useful 
in perturbative calculation, loses its position because the dynamical symmetry breaking 
appears as a nonperturbative phenomenon. Nevertheless if we can use it even for non- 
perturbative calculation, we may expect that computation will become easier and the 
relation between perturbative and nonperturbative phenomena will become clearer. In 
this paper we introduce the effective Langevin equation for two-dimensional nonlinear 
sigma model in the large N limit. With the help of it we derive some nonperturbative 
results which, of course, coincide with the results of the path-integral quantization. 

We consider the two-dimensional nonlinear sigma model whose Lagrangian in Eu- 
clidean space-time is 

L = T^dpfcdpfc (1) 



with the constraint 

N 

<j>i<fii = -o- (2) 

g 2 

Introducing a Lagrange multiplier field a and taking account of the constraint, we 
rewrite the Lagrangian as 



1 IN 

L c = A<f>i d »<i>i + o^Wi^i j)- ( 3 ) 



In the stochastic quantization scheme fields depend on the fictitious time t and 
their development along the fictitious time obeys the Langevin equation 



dL 

<t>i (x, t ) = - — ± I ^ xf ) +Vi(x,t), (4) 



where a dot denotes fictitious-time derivative and r)i(x,t) is a noise field whose corre- 
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lations are 



{rn(x,t)) ri = 0, (5.a) 

(r)i(x, t^x', t% = 2S ij S(x - x')S(t - t'). (5.6) 
Here the average manipulation (■ • -) n is defined as 

{"\ = \ J ^(•••) ex p{-^ j d 2 xdtrn(x,t)rii(x,t)y (6.a) 

Z = y _D?yexpj — - J d 2 xdti]i(x,t)rii(x,t)^ . (6.6) 

The Langevin equation for the field fa with the Lagrangian (3) becomes 

4>i(x, t) = d 2 4>i(x, t) - a(x, t)4>i(x, t) + r]i(x, t). (7) 



The equation (7) is studied perturbatively. [51 In that case the fields fluctuate around a 
point on the (N — l)-dimensional sphere (<f>i<f>i = jp), so that the vacuum loses the 
original O(N) symmetry and, as a result, (N — 1) independent 0's become massless 
Nambu-Goldstone bosons. Nevertheless these particles suffer from infrared divergence 
in two dimensional space time. On the other hand it is well-known that if we study 
the model nonperturbatively in the large iV limit by ordinary quantization methods, 
the O(N) symmetry recovers and the field gains a dynamical mass! 61 In the following 
we introduce the effective Langevin equation for the two-dimensional nonlinear sigma 
model in the large iV limit and investigate nonperturbative phenomena. The method 
used here is applicable to other models with dynamical symmetry breaking. 

To evaluate a(x, t) in the Langevin equation (7) we regard it as a functional of 
<f>i(x,t) and insert the unity 



/ 



D(f> | J | 5{<pi(x, t) - d 2 (f)i(x, t) + a(x, t)(f>i(x, t) - rji(x, t)} 
into the generating functional (6.b) and integrate over the noise field r\. Here | J 
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denotes the Jacobian factor which is a divergent constant ^ . We obtain 



Z = J D(j>\ J \ ex p{-^ J d 2 xdtLy (8.a) 



where 



L = hj>i(x, t)-d 2 <Pi(x, t) + a(x, t)<t)i{x, t))(<j>i(x, t)-d 2 <Pi(x, t)+a(x, t)<j>i(x, t)). (8.6) 

We thus derive a field equation for a(x,t) 

<9L 

— = i (0 l - ( 9 2 i + a0 l ) = O. (9) 
oa 

Restoring the noise fields as the integration valiables of the generating functional (8), 

rgi 

we can obtain the Langevin equation (7) again, in which a(x, t) is the solution of the 
field equation (9) 

a 2 

a(x,t) = -^</>i(x,t)(<f>i(x,t) - d 2 Ux,t)), (10) 

where the constraint (2) is taken into account. In large iV expansion, since 5u = N, 
we may expand a(x, t) as 

a(x,t) = (a) v (t) + 0(l/N) + ---. (11) 

To observe dependence of {a) v on the stochastic time, we call the Langevin equation 
for it; 

. dL c 

N 

= - ((pi(x, t)(pi(x, t) - -=) + 

Noting that the expectation value of one noise should vanish and taking account of the 
constraint (2), the expectation value of the above Langevin equation vanish 

(c^C = o, 

which holds at any fictitious time. Consequently (a) v does not depend on the fictitious 
time and we will evaluate it at t — > oo later. In the following we study the Langevin 
equation (7) exclusively in the large N limit N — > oo, which enables us to treat it 
nonperturbatively since the Langevin equation (7) no longer has any interaction terms. 



f We use Ito's calculation rule in this paper. 
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To investigate the dynamical breaking of the O(N) symmetry, we decompose 
(f>i(x,t) into the expectation value {<f>(x, t)) v = $i(t), which is the order parameter 
of the O(N) symmetry, and a fluctuating field </>•(£, t). Then the Langevin equation 
(7) becomes 

$(x,f) = (d 2 - {a)M{x,t) - {$i(t) + (a) v ^(t)} + Vi (x,t). (12) 

(a) v plays a role of square of dynamical mass in the equation (12). By definition of 
(ft'^k = 0,t), its expectation value must vanish and consequently 

$i(t) = -<<7)„$i(f), (13) 

which explicitly shows a relation between the order parameter $j and the square of 
the dynamical mass (cr)^. If (a) v = 0, the equation (13) imposes no restrictions on 
Nevertheless it is easily known that the massless bosons cause infrared divergence. On 
the other hand if (a) v > 0, $j vanishes when t — > oo, namely, the O(N) symmetry 
does not break down. 

We easily obtain the solution of the equation (12); 

t 

ti(k,t) = J dre-^+^-^^r)}, (14) 

where 

#(M) = J ^e^ti&t), 

TH(k,t) = J ^e ikx Vi (x,t). 
The correlations of ^(/c,t)'s are 

(rn(k,t)) v = 0, (15.a) 

(7fc(M)»7#',0>»/ = 25 tj 5(k + k')5(t - t'). (15.6) 
Let us derive a gap equation and compute (0%. Using the relations (10) and (13), 
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(a) v is written as 



(0% =— lira (<f>i(x,t)(-fc(x,t) + <9 2 0j(x, t))) v 

1\ t— >oo 

=^ lim (ct^Xa)^) + ^(x,09 2 ^(x,0 + MS 2 (0))„. 

1\ t^oo 



(16) 



In the second equality we have taken account of ((//(&;, t))^ = and Ito's calculation 
rule 

»=fi»«*=a{s*+8& (0) }- 

We can easily compute the last term with the equations (14) and (15) 



^ lim (ti(x,t)d 2 ti(x,t)) v 

iV t— >oo 

g2 r aI^aIV 



I d -^ r e-^ k > lim (-MX, t)^(k, t)) v 



N „ 

■) i d k k 



= ~9 



• [Jo 

J (2^ 



(2tt)2 A;2 + (a), 

# 2 A 2 # 2 A 2 

(cr) v ln(— — ) 



4tt n v (aV 4tt ' 

where we have introduced a straight cutoff parameter A. The two divergent terms in 
the above equation and in the equation (16) cancel out. Interested in non-zero solutions 
of (a) v we divide the both sides of the equation (16) by it and obtain the gap-equation 



Note that the above equation can be obtained by computing the expectation value of 
N/g 2 = (&(x,t)&(a;,*)>„- 



| The equation (16) can be directly derived from the equation (7). 
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If we define the renormalized coupling constant g r at a renormalization point /i as 

ill, r A* 

the solution ot the gap-equation (17) is 

/ \ 2 r 47T 4^(00)^(00) 

(<r)„ = /x exp{--^ + }, (18) 

9r 

which is independent of the renormalization point fi. Remembering the relation (13), 
we conclude that the O(N) symmetry is not broken, that is, 

$i(oo) =0 

and the n scalar particles gain a mass m 

m 2 = /! 2 exp{--^}. 

9r 

These results are same as the path integral approach gives. 

In this paper we have not discussed stability of the vacuum. Nevertheless the 
methods developed here can be applied to, for example, the Nambu-Jona-Lasinio 

ryi 

model whose symmetry spontaneously breaks down in the large N limit. In that case 
the Langevin equations for bound states are introduced and stability of the vacuum is 
considered convergency of them. 
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